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Abstract: In this paper all eight-vertex type solutions of the col- 
ored Yang-Baxter equation dependent on spectral as well as color 
parameter are given. It is proved that they are composed of three 
t~^ ' groups of basic solutions, three groups of their degenerate forms 

Q,. ' and two groups of trivial solutions up to five solution transforma- 

'""' ' tions. Moreover, all non-trivial solutions can be classified into two 

7? ' types called Baxter type and Free-Fermion type. 

Q: 

lO ! §0. Introduction 

>: 

^D , The Yang-Baxter or triangle equation which first appeared in Refs. 1-3 plays a promi- 

^^ , nent role in many branches of physics, for instance, in factorized 5-matrices [4], exactly 

^^ . solvable models of statistical physics [5], complete integrable quantum and classical sys- 

f"^ . tems [6], quantum groups [7], conformal field theory and link invariants [8-11], to name 

-■««« . just a few. In view of the importance of the Yang-Baxter equation, much attention has 

>: 



% 



been directed to the search for solutions of the equation. 

Colored Yang-Baxter equation dependent on spectral as well as color parameters is 
a generalization of usual Yang-Baxter equation. It has also attracted a lot of research 
interest (see Refs. 12-16) to find exact solutions for this type of Yang-Baxter equation. 
This is because the colored Yang-Baxter equation concerns free Fermion model in magnetic 
field, multi- variable invariants of links and representations of quantum algebras and so on 
(see Refs. 17-21). 

The eight-vertex type solution of the colored Yang-Baxter equation has been investi- 
gated previously in Refs. 12, 17 and 20. In Ref. 17 Fan and Wu first provided a single 
relation between the eight vertex weights in the general eight-vertex model by the pfaffian 
or dimer method, the so-called free-fermion condition. Based on this work, V.V Bazhanov 
and Yu. G. Stroganov obtained a eight-vertex solution for the colored Yang-Baxter equa- 
tion in Ref. 12, devoted to the eight-vertex free fermion model on a plane lattice. In Ref. 
20 J. Murakami gave another eight-vertex solution in discussing multi- variable invariants 
of links. These are only two eight-vertex solutions for the colored Yang-Baxter equation 
we have known up to now. 
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The main theme of this paper is to give and classify ah eight-vertex solutions of the 
colored Yang-Baxter equation. The way to find the solution is from a computer algebra 
method given by Wu in Ref. 22. Moreover, a theorem in the Ref. 22 can prove that all 
solutions can be obtained by Wu's method. The paper is organized as follows. In section 
one we will review the colored Yang-Baxter equation which in fact is a matrix equation. If 
the equation is expressed in components form it is composed of 28 polynomial equations 
in the eight-vertex case. In this section we will first introduce the symmetries or solution 
transformation for this system of equations and some definitions of Hamiltonian coefficient, 
initial value condition, unitary condition and non-trivial gauge solution of the colored 
Yang-Baxter equation. Using the symmetries we can simply the system of equations as 12 
polynomial equations. In section two we will apply the computer algebra method to the 12 
polynomial equations to get the algebraic curves and differential equations satisfied by the 
eight-vertex type solution of the colored Yang-Baxter equation. In this section we will also 
give two relations satisfied by Hamiltonian coefficient which will play an important role 
in classification of eight-vertex solutions of the colored Yang-Baxter equation. Based on 
the second section, in the third section we will construct all non-trivial gauge eight-vertex 
type solutions of the colored Yang-Baxter equation and classify them into two types called 
Baxter and Free-Fermion type. The fourth section is devoted to general solutions and the 
relation between Hamiltonian coefficients and spin-chain Hamiltonian. In the third and 
fourth sections we will also show the two solutions appeared in Refs. 12 and 20 are special 
cases of the general solutions obtained in this paper. 

In this paper symbolic computation will be applied to accomplish some tedious com- 
putations and the results obtained by computer calculations will be denoted by the symbol 



§1. Colored Yang-Baxter equation, its symmetry and initial 
condition 

By colored Yang-Baxter equation we mean the following matrix equation 

Ri2{u, C, r])R23{u + V, C, X)Ri2{v, 7], A) = R23{v, V, A)i?i2(n + V, C, X)R23{u, C, v), 

Ruiu, S„ v) = Riu, C, v) «) E, R23iu, S„v) =E(S) Riu, C, i]), 

where R{u, £,, if) is a matrix function of A^^-dimension of u, $^ and r/, E is the unit matrix of 
order N and dS) means tensor product of two matrices, u, v are called spectral parameters 
and ^, rj colored parameters. If the matrix is independent of colored parameters, then 



the colored Yang-Baxter equation (1.1) will become as the usual Yang-Baxter equation. 
If it is independent of spectral parameter then (1.1) will be reduced to the pure colored 
Yang-Baxter equation. 



Rui^, V)R23{t ^)Rl2iV, A) = R23ir], X)Rl2iC, A)i?23(^, V)- 



(1.2) 



For the colored Yang-Baxter equation (1.1), the main interest in the paper is to discuss 
the solutions with the following form 

( Rll{u,C,r]) Rll{u,C,v) \ 

Rlliu,C,v) Rhiiu,C,v) 

RlUu,C,v) R2iiu,C,v) 

\ Rll{u,tv) R^Uu^tv) J 



R{u,^,r]) 



:i.3) 



The eight weight functions in (1.3) are denoted by 

ai{u,tr]) = Rll{u,^,'n), 05(^,^,1?) = Rll{u,^,r]), 
a2{u, t V) = Riliu, C, V), aeiu, ^, r?) = Rl^{u, ^, i]), 
a3{u,tv) = R2i{u,^,v), ajiu,^,!]) = R\l{u,i,r]), 

ai(u, ^, f]) = R22{u, ?, v), as{u, t v) = Rii{u, t v)- 

We call the solution (1.3) eight-vertex type solution if it satisfies in addition ai{u, S,,7]) ^0 
(i = 1, 2, • • • , 8). Further, we only consider the solutions which are meromorphic functions 
of u and ^, rj. For notational simplicity, throughout this paper let 



Ui = ai{u,^,r]), Vi = ai{v,r],X), Wi = ai{u + v,^,X), 



l,2,---,8. 



For eight-vertex type solutions, the matrix equation (1.1) is equivalent to the following 28 
equations: 



urw^vs - usW2Vy = 0, 

U7WSV3 - UsWyV2 = 0, 

U2WSV2 - U3W2V3 = 0, 
U2WSV7 - u^wjvs = 0, 



(1.4a) 



U1W5V2 + ujwsvq — V2W1U5 — V5W2U3 = 0, 

UiWlVj + UJW3V4 — VjW^U^ — VlWjU^ = 0, 

U2WQV1 + U5W7VS - VQW1U2 - V3W2UQ = 0, 

U1W2V1 + U7W4V8 — V2W1U2 — V^W2U% = 0, 
UiWjV^ + UjWqVs — VTW^U2 — ViWjUq = 0, 
Ui'WjV2 + UjWqVq — ViWiUi — ^7^2^*4 = 0, 

U4WQV2 + ujwsv^ - V2W4UQ - VQW2U3 = 0, 

t(4W4f 7 + U-jW^Vl — VjWqUq — VaWiUz = 0, 
U2'W5V4 + UqWjVs - V5W4U2 - V3W2U5 = 0, 

U4W2V4 + urwivs - V2W4U2 - VQW2U5 = 0, 

UiWjVQ + U7W5V3 — VjWqU2 — V4W7U5 = 0, 
UiW^V2 + UJW5V5 — V4W4U7 — VjW2Ui = 0, 

U1W5V3 + UsWjVq - V3W1U5 - V5W3U2 = 0, 

uiwivs + USW2V4 - VSW5U5 - V1WSU2 = 0, 
U3WQV1 + u^wsvr - VQW1U3 - V2W3UQ = 0, 

U1W3V1 + USW4VJ — V3W1U3 — V^WsUq = 0, 

U1W8V5 + U8WQV2 - VSW5U3 - viWsUQ = 0, 

U1WSV3 + UsWqVq - ViWiUs - VSW3U4 = 0, 



(1.4b) 



(1.4c) 



(1.4d) 



(1.4e) 



U4WQV3 + USW7V5 — V3W4UQ — VQW3U2 = 0, 
U4W4VS + USW2V1 - vsvjQUe - V4WSU2 = 0, 

U3W5V4 + UqWsV7 - V5W4U3 - V2W3U5 = 0, 
U4W3V4 + U^WlVj — V3W4U3 — VqW^U^ = 0, 

U4W8V6 + USW5V2 - vsweU3 - V4W8U5 = 0, 

U4WSV3 + USW5V5 — V4'W4Us — VSW3U1 = 0. 

Assume R{u,^,rj) is a solution of (1.1). Having carefully studied the system of equa- 
tions (1.4), we find there are five symmetries for eight-vertex type solutions of the colored 
Yang-Baxter equation (1.1). 

(A) Symmetry of interchanging indices. 

The system of equations (1.4) is invariant if we interchange the two sub-indices 2 and 
3 as well as the two sub-indices 7 and 8 or the sub-indices 1 and 4 as well as the two 
sub-indices 5 and 6. 

(B) The scahng symmetry. 

Multiplication of the solution R{u,^,ri) by an arbitrary function g{u,S,,rj) is still a 
solution of the colored Yang-Baxter equation (1.1). 

(C) Symmetry of weight functions 

If the weight functions a2{u,S,,7]),as{u,S,,T]),aj{u,S,,7]) and as{u,^,ri) are replaced by 
the new weight functions 

respectively or arj{u,^,ri) and aQ{u,^,ri) by —a5{u,^,7]) and —aQ{u,^,7]), where N{^) is 
an arbitrary function of colored parameter and s is a complex constant, the new matrix 
R{u,S,,rj) is still a solution of (1.1). 

(D) Symmetry of spectral parameter. 

If we take a new spectral parameter u = fiu where /i is a complex constant independent 
of spectral and colored parameters , the new matrix R{u',^,rj) is still a solution of (1.1). 

(E) Symmetry of color parameters. If we take new colored parameters C = fiOi 
9 = f{ri), where /(^) is an arbitrary function , then the new matrix R{u,C,0) is also a 
solution of (1.1). 



The five symmetries (A)-(E) are called solution transformations A-E of eight-vertex 
type solutions of the colored Yang-Baxter equation (1.1) respectively. 

Dividing both sides of the third of (1.4a) by a2(n, S,, rj)a2{u + v, $,, X)a2{v, rj, A), we get 

f{u + v,tX)=f{u,^,v)f{v,V,>^), (1-5) 

where f{u,^,7]) = a^i^u, £,, rj) / a2{u, ^, rj) . Putting u = v = r] = in (1.5) we have 

/(0,e,A) = /(0,e,0)/(0,0,A). 
Substituting this formula into the one obtained by taking u = v = ^ = in (1.5) we get 
/(0, 0, A) = /(O, 0, ry)/(0, r?, A) = /(0, 0, r/)/(0, r/, 0)/(0, 0, A). 

This means 

/(0,r/,0)/(0,0,r/)=l. 

Otherwise, it is easy to show that f{u,S,,rj) = 0, i. e. a^^u,^^,!]) = 0. Therefore 

/(0.e.,) = ||, (1.6) 

where M(^) = /(O, ^, 0). On the other hand, if we differentiate both sides of (1.5) with 
respect to the spectral variable v and then set v = 0, X = r], then 

/'(u,e,r/)=/'(0,r/,r/)/(n,^,r?) (1.7a) 

holds, where the dot means derivative with respect to u and the simple formula 

dG{u + v) _ dG{u) 
Tv '^=0" du ' 

for any function G{u), is used. Similarly, one also has 

/>,e,A) = /'(0,e,O/(«,e,A) (1.7&) 

if we differentiate both sides of (1.5) with respect to u and then set u = and i] = ^. 
The two formulas above imply f'{0, £,,£,) is a constant independent of colored parameter 

^. Hence 

M(£) 
f{u,C,v) = j^eM^u), (1.8) 

where z^ is a complex constant. 



From the group of equations (1.4a) we have 

f{u, ^, 'n)h{v, 7], A) = h{u + V, ^, A), (1.9a) 

f{u, e, v)f{v, V, A) = f{u + V, C, A), (1.9b) 

h{u, C, v) = Hu + V, ^, X)f{v, r], A), (1.9c) 

Hu, e, V) = fiu + V, e, X)hiv, f], A), (1.9d) 

where f{u,(,r]) = as{u,^,r])/a2{u,(,r]), h{u,^,r]) = as{u,^,r])/a7{u,(,r]). If we let v = 
and A = in (1.9d) then 

/i(n,e,r/)=/(u,^,0)/i(0,r?,0). 

Substituting this into (1.9c) and using (1.9b), one obtains 

f{u, 1 0)h{0, 7], 0) = f{u, C, 0)f{v, 0, 0)/i(0, A, 0)f{v, r?, A) 

or 

h{0,v,0) = f{v,0,0)h{0,X,0)f{v,7],X). 

This formula imphes u = in (1.8) and then one can obtain 

where I is a complex constant independent of spectral and colored parameters. 
So, up to the solution transformation B and C one can assume 

asiu, C, v) = a2{u, ?, rj) = 1, as{u, ^, r?) = 07(14, ^, r?) 

without losing generality. In the case (1.4a), (1.4b), (1.4c) and (1.4d) can be simplified to 
the following 12 equations. 



(1.10) 



a5{v,ri,X) + a5iu,^,ri)ai{u + v,^,X) 

-ai{u,(,r])a5{u + v,^,X) - a7{u,i,r])aj{u + v,^,X)aQ{v,'q,\) = 0, 

a7(n + V, ^, \)ai{v, t], A) + 05(1*, ^, r])a5{u + v, ^, X)a7{v, t], A) 
-ai(u,^,r/)ai(n + t',^,A)a7(t;,r?,A) - a7{u,^,r])a4{v,7], X) = 0, 

a6(w,^,??) +ai(u + f,^,A)a6(i',??,A) 

-aeiu + v,^, X)ai{v,r], X) - a5{u,£,,r])a7{u + v,C, X)a7{v,r], X) = 0, 

a6{u,^,r])a5{v,r],X) +ai{u + v,^,X) 

-ai{u,(,r])ai{v,r],X) - a'r{u,C,r])a4{u + v,^,X)a7{v,r], X) = 0, 

aeiu, C, r])a7{u + v, C, X)ai{v, rj, A) + a^iu + v, ^, A)a7(w, r?, A) 
-ai{u,(,r])a7{u + v,^,X)a5{v,r],X) - a7{u,£,,r])a6{u + v,C, X) = 0, 

aviu, t V)ai {u + v, C, X)ai {v, rj, A) + 04 (u, (, v)a7 {v, rj, A) 

-ai{u, C, 'n)a7{u + w, ^, A) - 07(1^, ^, r?)a6(u + v, C, X)aQ{v, r/, A) = 

plus six equations obtained by interchanging the sub-indices 1 and 4 as well as 5 and 6 in 
each of Equations (1.10). We call the six equations counterparts of (1.10). 

Now we solve the equations obtained by letting u = and ry = ^ in (1.10) with respect 
to the variables {oi(0,^,^),a4(0, ^, ^), a5(0,^,^),a6(0, ^, ^), 07(0,^,^)}. It is easy to prove 
that 

Proposition 1.1 For a solution of equations (1.10), weight functions satisfy 

ai(o,e,e) = 04(0,^,0 = 1, 

(1.11) 

«5(o, e, = «6(o, e, = "7(0, e, = 08(0, e, = 0. 

Otherwise, up to the solution transformations A, B, C, D, and E, we have two trivial 
solutions of the colored Yang-Baxter equation (1.1). The first is 

(1.12a) 
a2{u,tv) = a3{u,tr]) = ajiu,^,-!]) = as{u,^,r]) = 1, 

where H{u, S,, rj) is an arbitrary function of spectral parameter u and colored parameters 
^, r] and the second 

F(0 
ai{u,i,r])=aA{u,i,ri)=arj{u,i,'n) = -a(i{u,^,r]) = ——-exp{u), 

F{v) {1.12b) 

ai(.u, ^, f]) = a2{u, ^, f]) = 1, ariu, C, v) = M^, ?, v) = h 



where i? = —\ and -F(^) is an arbitrary function of colored parameter ,^. 

Definition 1.2 By a gauge solution of the colored Yang-Baxter equation (1.1) we 
mean the solution whose weight functions satisfy a2{u,^,7]) = a^{u,^^r]) = 1 and 
ar{u,^,ri) = as{u,^,r]) and the condition (1.11). 

(1.11) is called an initial condition of gauge solutions. The condition is simple but 
very important. It will be quoted again and again in finding gauge solutions of the colored 
Yang-Baxter equation. For example, taking v = —u, A = ^ in (1.10) and their counterparts 
and using the initial condition (1.11), one has the following equations. 



(1.13) 



a6{-u,r],0 = -a(i{u,^,ri), 

a5{-u,ri,^) = -a5(u,C,^), 

1 -ai{u,^,r])ai{-u,r],C) + a6{u,(,r])a5{-u,ri,^) - 07(^,^,7^)07 (-n, 7?,^) = 0, 

1 -a4(n,^,r?)a4(-n,77,^) + a5{u,(,r])aQ{-u,r],^) - a7(u,^,r/)a7(-n, 77,^) = 0. 
The unitary condition of a solution of Yang-Baxter equation means 

i?(0, e, = E, Riu, C, f]) Ri-u, f], = 9iu, C, V) E, 

where E is the unit matrix and g{u, ^, rj) a scalar function. Hence, it is easy to get from 
(1.13) that 

Proposition 1.3: For gauge solutions R{u, (,, rj) of the colored Yang-Baxter equation 
(1.1), the unitary condition is 

R{u, ^, 7]) R{-u, 7], ^) = (1 - a^iu, ^, 7]) aeiu, ^, r?)) E. 

Differentiating both sides of all equations in (1.10) and their counterparts with respect 
to the variable v and letting f = 0, A = 77, by virtue of the initial condition (1.11) one 



immediately obtains 

m^ir]) + a5{u,^,r])a[{u,^,r]) - ai{u,(,r]) a'^{u,^,r]) - mQ{r]) a7{u,C,r])'^ = 0, 

a'7{u,i,r]) + {mi{r]) - mi{ri))a'j{u,i,r]) + 7717(7]) {a5{u,S,,vf - ai(.u,S,,vf) =0, 

a'siu,^,?]) -mQ{r])ai{u,^,r]) + mi{r]) aQ{u,^,r]) + 171,7(77) 05 (n,^,r?) 07(^,^,7/) = 0, 

a'i{u,^,r]) -771,1 (77) ai(n,^,r?) + ^5(77) a6(n,^,r7) - 771,7(77) a4(n,^, 17) 07(^,^,7/) = 0, 

ae{u, ^, 77) a'jiu, C, v) - aiiu, ^, 7/) Qq^u, ^, 77) 

+7111(77) a6(n,^, 77) 07 (u,C,^) -7ii5(r/)ai(7x,^,7/)a7(n,^,77) +7717(77)05(14,^,77) = 0, 

a7(u, C, 77) a'l (71, ^, rj) - ai {u, ^, 7/) 07(1*, ^, 77) 

+mi(r/)ai(n,^, 77) 07(7^,^,7/) - meir]) a6{u,(,r]) a7{u,^,r]) +7717(77)04(14,^,77) = 

(1.14a) 
and their counterparts, where and throughout the paper we denote 

d 

for i = 1,4,5,6,7. 

We call mi{^) Hamiltonian coefficients of weight functions with respect to spectral 
parameter or simply coefficients. Sometimes we write 777j instead of mi{^) for brevity. 

If we differentiate (1-10) with respect to u and let 71 = 0, 7/ = ^ and then replace the 
variables v and A by n and r/, we have 

m^si^) + a(,{u,i,vi) a']^{u,£,,r]) - ai(7x,C,i/) ^(tx,^!/) - m^{i)a-j{u,^,r]f = 0, 

a'liu, C, V) + ("11(0 - iTT-iiO) ariu, ^, 77) + 7717(0 {aeiu, C, v)"^ - ai{u, (, rjf) = 0, 

4(^>C,^) -"^5(001(^,^,77) +777,1(005(74,^,77) +777,7(0 06 (7i,0??)«7(^^,0^) = 0> 

a'i{u,i,rj) -r77i(Oai(n,^,77) + 777,6(0 05(^^,0??) - "1-7(0 a4(«,0??) 07(^,0??) = 0> 

a^{u, ?, r]) 07(77, ^, rf) - aj{u, ^, 77) 05(71, ^, 77) 

+mi{i)a^{u,i,ri)a7{u,i,ri) - nieiO ai{u,(,r]) a7{u,(,r]) + "^7(0 a6(w,C,^) = 0, 

07(77, ^, 7/) a'i(77, ^, 7/) - ai(7i, ^, 77) o'7(7i, S,, 77) 

+mi(0ai(ii,?,i/)a7(u,?,^) - m^i^) a5{u,(,r]) a7{u,(,r]) + ?t77(0 «4(ii,C,^) = 

(1.146) 

and their counterparts. 
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Remark 1 : In this paper , the following trick will often be employed to obtain the 
equations like (1.14a) and (1.14b): We first do the calculation with respect to v and take 
A = ?7 to yield some equations. Then repeat the same operation with u in place of v and 
take ?7 = .^ to yield another equation. Then we compare the two results to reduce to some 
of the formulas. For example, formula (1.8) is obtained by this method. 

Comparing (1.14a) with (1.14b), we found that the trick, in fact, is to interchange the 
subindexes 5 and 6 (or 1 and 4 ) and then replace mi{ri) by mi{^) in original equations 
(1.14a). We call this trick symmetric operation . 

Remark 2: If the variable with respect to which we differentiate (1.10) is not a 
spectral parameter v but a colored parameter A, the equations obtained are the same as 
(1.14a) if we still let f = and X = rj. Of course, then the dot means the derivative with 
respect to the colored parameter ry, that is the second colored variable in ai{u,^,r]), and 
mi{rj) = dai{v,r],X)/dX\^^Q^ . Similarly, (1.14b) also represents the equations obtained 
by differentiating (1.10) with respect to the first colored variable S, in ai{u, ^, rj), but in this 
case the dot means the derivative with respect to ^ and mi{^) = dai{u, ^, 7]) / dX\^^Q c- 

§2. The coefficients, curves and differential equations of weight| 
functions 

In this section we will discuss properties of Hamiltonian coefficients, curves and differ- 
ential equations satisfied by weight functions. 

It follows from the second equation of (1.14a) and its counterpart that 

2a7(u,C,^) + 1717(1]) {a5{u,^,r]f +aQ{u,^,r]f -ai{u,S,,vf -a4{u,S,,vf) = 0. (2.1) 

If the symmetric operation is used then one also has 

24(u,e,r/) +m7(0 {a5{u,C,r]f + ae{u,C,r]f - ai{u,C,r]f - aA{u,C,r]f) = 0. (2.2) 

We know (lig + Ug — uf — U4) ^0 due to the initial condition (1.11). Compared the formulas 
(2.1) and (2.2), one shows 

mriO = mj{r]). 

If rm{C) = then 

ariu, 1 1]) = 0, mi (77) = 7714(77) 

due to the second of (1.14a) and its counterpart. This implies ai{u,S,,i]) is a function of 
only the colored variables ^ and rj. 
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Proposition 2.1: For gauge eight-vertex solutions, mi{^) is a constant independent 
of colored parameters and aj{u,S,,rj) is independent of spectral parameter if mi{ri) = 0. 

Remark 3: In fact, as mentioned in remark 2, the latter property of proposition 
2.1 also holds for colored parameter, namely, aj{u,S,,ri) will be independent of the colored 
parameter ^ (or ij) H d/d^{a7{u,^,r]))\^^Q^^^i: = (or d/dr?(a7(u,C,^))L=o,c=r, = 0)- 

In what follows we denote mi{S,) = a. 

Furthermore, if a = 0, letting u = 0, ry = ^ in the following equation 

o'i{u,i,r]) =m6{r])a6{u,^,r]) - mi{r])ai{u,^,r]), (2.3) 

which is from the sixth in (1.14a), we obtain mi{rj) = owing to the initial condition 
(1.11). Therefore m4(^) = 0. Similarly, we have 

«6("i^'^) =Tm,i{r])a(i{u,i,r]) - m^{r])ai{u,i,r]) (2.4) 

from the fifth in (1.14a) and then 

"^6(??) = -m^iv)- (2-5) 

Substituting 171^(7]) = —m^{rf) and mi{ri) = m^^r]) = into the third and fourth of (1.14a) 
and their counterparts, we have the following differential equations. 



Therefore, if a = then weight functions satisfy 



(2.6) 



—^ai{u,^,ri) =1715(7])'^ ai{u,^,ri), z = 1,4,5,6. (2-7) 

Furthermore, using the symmetric operation 

—^ai{u,^,7i)=77i5{^fai{u,^,7i), i = 1,4,5,6 (2-8) 

hold. Hence 771,5(77) actually is a constant independent of colored parameters and not 
identically zero, otherwise, the solutions will be independent of spectral parameter. Thus 
we can let 771.5(77) = (3. 

The argument above implies the following proposition. 

Proposition 2.2 For a gauge eight-vertex solution, there exists at least one between 
77ij{^) and 7775 (^) (or 77ig(^)) which is not zero identically. Otherwise, the solution will be 
independent of spectral parameter. 
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As for the Hamiltonian coefficients mi{^) and w,4(^) we have 

Proposition 2.3* : For gauge eight- vertex solutions of the colored Yang-Baxter equa- 
tion (1.1) 

miiCf - m^iCf = 0. 

Proof: As first step, we regard the weight functions Wi {i = 1, 4, 5, 6, 7) in (1.10) and 
their counterparts as indeterminates, The left side of each of (1.10) and its counterpart 
are polynomial functions of the indeterminates. After eliminating the five indeterminates 
{wi^Wi, W5,wq, itJy} in numerically increasing order with respect to the system of equations 
(1.10), we can obtain seven equations which do not contain the indeterminates Wi {i = 
1,4,5,6,7). Then differentiating them with respect to the spectral variable v , letting 
u = 0, A = r/ and then substituting the initial values (1.11) into the resulting ones, we 
obtain the following seven polynomial equations 

myiif — nijUiul — 3miUiU7 + m^^uiuj — m-ju^Uj — ni'jUi + m^UQUj + tuqUqUj = 0, 

—niQUiu^ + nijUiUQU-j + ttiyu^u^uy + miU4UQ + m^^u^UQ — itiqu^uq — m^Uj + m^ = 0, 

m-ju\ — m'ju\ — m-jul + itijUq + 2 m^^uj — 2 niiuj = 0, 

—m^uiu^ + niiuiu^ + m^uiu^ + m-juiu^u-j + m-ju^^u^uj — m^u^UQ — m^Uj + m^ = 0, 

mjUiUQ — mQUiWj — m^uiui — miulu^ + mju^Uj + m-ju^ + m^UQWr + miUQUj = 0, 

m-ju\u^ — rriQUiU/^wj — m-juiu^ + 2 7714^11x51x7 — 2miUiu^U'j + mjUiUQ 
—mjU^u^Uj + m^u^UQUj + itlqUj — m^uj = 0, 

—m-iUiUi + nijUiUj + nijui + mju^ul + m^UiUj + miu^ui — m^u^uj — uiqu^uj = 0. 

(2.9) 
As second step, we think of rrii as indeterminates and first eliminate mi, 7714 and 777.5 to 
get two systems of equations, which are equivalent to (2.9). The first is 

— 77757X1U4 -|- 7771^17X5 -|- 77747717x5 + m-jUiUQU-i + myu^u^uj — m^u^UQ — uiqu^ + 7775 = 0, 

—m-ju\u^ + m-jUiulu^ + Im^uiu^uj + m-jUiu^ — nijUiu^UQ — Am^uiu^uj 
—2 mjUiUQUj — 2 nijUiU^Uj + 2 m^u^u^u-j + 2 m^u^ — 2 77757x7 = 0, 

—7)777(17x47*5 + meuiu^uj + mjUiu^UQ — m-jUiUQ + m-jU/^u^v^ — m^u^UQWj 
—niQUj + 777,57(7 = 0. 

(2.10) 
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(2.11) 



which contains rrii , m4 and m^ . The second is 

{uiU4 + u^uq — Uj — l) {—m-juiu\u^ + rriQUiUiUi + m-juiu^u1 
—mjUiUQ + m-ju^u^ — m^u^u^uj) = 0, 

ui (uin4 + u^uq — Uj — 1) {mjuiu^ul — m^uiu^u-j — mju\u^UQ + niQUiUQUj 
-mjuluQ + m,ul + m^u^ul - m^ul) = 0, 

ui [uiUi + tisue — U7 — l) (— m7«4«5t(6 + m%u\u%ui + m7n4ti5ti7 + mju^u^UQ 
—mju^UQ — mQuluQUj — m^u^u^ + ttiqu^u-j) = 0, 

ui (tiin4 + tisue — U7 — 1) (—mjUiu^UQ + m-juiu^ — m-ju\u^ + ijiquIuj 
+mjU4ul + mjUiU^UQ — mju^UQ — niQuluj) = 0, 

which do not contain mi,m4 and mr,. So the free fermion condition[17] 

U1U4 + u^uq — 1 — Uy = (2-12) 

or 

—mjuiu\u^ + m^uiUiUj + mjUiu^UQ — m^uiUQ + mju^u^ — niQU^ueUj = 0, 



+miU5UQ — mjUQ = 0, 



(2.13) 



—mQU^Uj + niQU^uj = 0, 

—rrijUiuluQ + m-jUiu^ — nijufu^ + 777,6^4^7 + m-jU/^u^ + m7ti4n5ng — m-jU/^UQ 
—m^u'lu'j = 

will hold. In the third step, applying the fourth in (2.13) as a main equation to kill the 
indeterminate rriQ in the three ones remained in (2.13) and then doing factorization of the 
new polynomial equations after killing Trig, one can obtain that 

—m'rUiuluQ + m-juiu^ — m,yu\u^ + m%u\u-j + mjUiU^ + mju^u^UQ — mjUiUQ 
-m6ulu7 = 0, 

mj {u^uq — 1) (nf ti5 — 2 uiu^uq + n^us — Ug + u^Uq) u^ = 0, 

(2.14) 
mj {U5UQ - 1) • 
{—uiuIuq + uiu^ue + uiu^Uj — uiu^ + n^us — 7*4^5 — u^uqUj + u^uq) u^ = 0, 

m-j {u^uq — 1) (— tifn4 + 2tiiti5n6 + uf — u^ul — u^Uq) us = 
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is equivalent to (2.13). It follows from the first of (2.14) that rrif ^0. Otherwise, thanks 
to proposition 2.2 04(^,^,7/)^ = a^{u,S,,ri)'^- The latter is impossible thanks to the initial 
condition (1.11). Hence the following three equations 

ufu^ — 2 U1U4UQ + u^tis — U5 + u^Uq = 0, 

— UiuIuq + UiU^Uq + UiU^Uj — U1U5 + tifus — U^ul — U^UqUj + U4UQ = 0, (2.15) 

— tif U4 + 2 U1U5UQ + uf — Uiul — U4UQ = 

and the first of (2.14) is equivalent to (2.14), where we use the initial condition (1.11) 
again to yield 1 — u^uq ^0. 

Finally, if we differentiate (2.12) and the third equation of (2.15) with respect to u 
and let u = 0, .^ = ?7 and then apply the initial conditions (1.11) again , we come to the 
conclusion of proposition 2.3. 

Remark 4: When we do the operation of eliminating indeterminates with respect to 
a system of polynomial equations, according to the theorem of zero structure of algebraic 
varieties [22], the coefficient of the term with the highest degree of the indeterminate in 
main polynomial equation (to be eliminated in other polynomials) should not be identified 
with zero. In the event it is identified with zero, we should add the coefficient into the 
system of equations to produce a new system of equations. Otherwise, it is possible to 
lose some solutions. For example, when we use the fourth equation in (2.13) to eliminate 
rriQ in the three remaining ones in (2.13), because the coefficient of mg in the fourth one of 

(2.13) is ^7(^4 — n^) which does not identify with zero due to the initial condition (1.11), 

(2.14) is equivalent to the system of equations (2.13). 

From the argument of proving proposition 2.3 above we see the system of equations 
(2.9) is equivalent to two groups of equations. The first is (2.10), (2.15) plus the first of 
(2.14). The second is (2.10) and (2.12), the free fermion condition. 

Now we consider the two cases respectively. For the first case we differentiate the 
second equation in (2.15) and take u = 0, ^ = rj. Then we substitute the initial condition 
(1.11) into the result to get 

mrj{ri) = m6(r/). (2.16) 

By doing factorization of the equation obtained by eliminating ii4 in the third equation of 

(2.15) and by using the second equation of (2.15) , we can get 

2 Uq (uq - U5) (uq + U5) {Ui - U5) {Ui + U5) {ui - Uq) {ui + Uq) = 0. (2.17) 

Together, (2.17) and the initial condition (1.11) imply 

aniu, C, v) = aeiu, ^, r?) (2.18a) 
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or 

abiu, t V) = -aeiu, ?, v) {2.18b) 

hold. Substituting (2.18a) into the third equation in (2.15), we then see 

ai{u,(,v) =a4{u,(,'n). (2.19) 

Substituting (2.18a) and (2.19) into the first equation in (2.14), we have 

(u5 — ui(u5 + ui){au5Ui — rriQUj) = 0. (2.20) 

So 

a ai{u, C, 7]) 05 (u, £,, 7]) - mQ{r]) ayiu, ^, 7]) = 0, (2.21) 

here the initial condition (1.11) is used again. (2.21) implies ra^{rj) ^0, or a will also 
identify with zero. But this will contradict proposition 2.2. If we do the symmetric 
operation with respect to (2.21) and use (2.19), then 

aai{u,^,rj)a^{u,i,ri) - mG{^) a7{u,(,r]) = 0, (2.22) 

(2.21) and (2.22) will imply niQ^rj) is also a constant independent of colored parameter. 
We let m^irj) = j3. 

If a^{u,S,,ri) = —aQ{u,S,,rj) we should have rriQ^rj) = and (2.21) still holds. Then 
a = 0. It is clearly impossible thanks to proposition 2.2. 

Combining (2.18a), (2.19), (2.21) with the third equation in (1.14a) we obtain 

{4{u, t V)? = /?' - (/5' - mi{i^f + a^) a^{u, ^, i^f + a^ a^iu, C, r,)\ (2.23) 

Using the symmetric operation we can show mi{r]) is also a constant independent of colored 
parameter. Let mi{rj) = 7. Similarly, 

{a[{u,i,il)f=l3^-{l3^--i'' + a^)a^{u,i,i^f + a^ai{u,i,r^Y (2.24) 

holds if we combine (2.18a), (2.19), (2.21) with the fourth in (1.14a). Substituting (2.16), 
(2.18a), (2.19) and (2.21) into the second of (2.10) we fund that the algebraic curve satisfied 
by the weight functions ai(ii, ,^,ry) and 05(^,^,7/) is 

a^ ujul -0^ul- (3'^ul + 2/37 U1U5 + /3^ = 0. (2.25) 

From the second group of equations, i.e. (2.10) and (2.12), it is easy to obtain 

l + a7{u,(,r]f -ai{u,(,r])a4{u,^,r]) - a^{u,i,ri)aQ{u,i,ri) = 0, 

a{ai{u,i,'q)a(i{u,i,ri) + a4{u,(,r])a5{u,^,r])) = [m^ir]) + meir])) a7{u,£,,r]), (2.26) 

a(ai(u,C,r/)2 +06(u,^,r/)2 - a4{u,(,r])'^ -05(^,^,7/)^) = 4mi(r?) 07(^,^,7/). 
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Then mi(rj) + 1714(1]) = 0. If we do the symmetric operation with respect to the third 
equation of (2.26) then we also obtain 

a{ai{u,C,r]f +a5{u,(,r]f - ai{u,i,'qf -aQ{u,(,r]f) = 4mi(^) 07(^,^,7/). (2.27) 

Letting ry = ^ in the third equation in (2.26) and (2.27) we reduce to 

ae{u,C,Cf = a5{u,(,0^- (2-28) 

Therefore, (2.16) and (2.28) give the following proposition. 

Proposition 2.4: For gauge eight- vertex solutions of the colored Yang-Baxter equa- 
tion (1.1) the Hamiltonian coefficients m^{ri) and niQ^rj) satisfy 

m5{i]f - mQ{rif = 0. 

From (2.26) and the second equation in (1.14a), we can calculate the weight function 
07(^,?,^) obeys 

7V is, /; 1^2.29) 

= a^ - {{m^{ri) + mfs{ri))'^ - 4ml{i]) - 2a^)aj{u,$,,rif + a^aj{u,$,,i]Y. 

Furthermore, as we said in remark 1, we can show {m^ijff' + mQi]))'^ — 2mi(ry)^ is a 
constant independent of colored parameter using the symmetric operation. Let S'^ = 
{m^{ri)'^ + mQ{r]))'^ - 2mi(r/)^. 

We conclude this section by the following theorem. 

THEOREM 2.5*: For a gauge eight-vertex type solution, its weight functions must 
satisfy one of two systems of equations. The first is composed of (2.18a), (2.19), (2.21), 
(2.23), (2.24) and (2.25). The second is composed of (2.26) and (2.29). 

§3. Gauge eight-vertex type solutions 

In this section sn{() and cd(C) = cn{()/dn{C) are Jacobian elliptic functions. 

Now we describe how to write down all gauge solutions of eight-vertex type of the 
colored Yang-Baxter (1.1) and classify them into two types called Baxter type and Free- 
Fermion type. 

3.1 Baxter type solutions. 

We consider the first case in theorem 2.6. Since the curve (2.25) only includes two 
weight functions ai{u,^,ri) and 05(^,^,7/) , we can parameterize ai{u,^,ri) and 05(^,^,7/) 
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as one-parameter functions. If /3 it a it 7 7^ 0, (2.23) and (2.24) are elliptic differential 
equations. Therefore, the solutions should be 

gn(An + F(0-F(r?)+/i) 
ai{u,i.,V) = a4{u,^,ri) = — — , 

a2{u,i,i) =az{u,i,r]) = 1, 

(3.1) 
, , . I , . , sn{\u + F{0-F{r^)) 

sn[fi) 
a7(n, ^, f]) = ±k sn{X u + F(^) - F(r/)) sn{\ u + F{i) - F(r]) + ^) 

where A;, as the modules of Jacobi elliptic function, is an arbitrary constant. 

If /3± a lb 7 = 0, the elliptic solutions (3.1) will degenerate into trigonometric solutions 

tan(A^ + F(g)-F(r/)+^) 
ai{u,i,,'n) = 04(^,^,77) = , 

a2{u,S,,v) = a3{u,S,,v) = 1> 

(3.2) 
. , ^ . , ^ , ton(An + F(0-F(r?)) 
a5{u,S.,v) = a(i{u,S.,v) = ± tonM ' 

a7(n,^,r?) = ±tan(An + F(^) - F(7?)) tan(An + F(^) - F(r/) +^). 

In (3.1) and (3.2) A ^ 0, /U 7^ are two arbitrary constants and F{S,) an arbitrary function. 

3.2 Pree-Fermion type solutions. 

Now we consider the second case in theorem (2.6). According to proposition 2.4 it is 
divided into two sub-cases, ra^{^) = mQ{S^) and m^{S,) = — w-6(C)- 

(3.2a) For the sub-case of m^iC) = w,6(^) 

Let m5(^) ^0 (we will put the case of m^{ri) = into the second sub-case ). It is 
clear that a ^ due to the second of (2.26). For brevity we let a = 1 up to the solution 
transformation D. When 

and 

the equation (2.29) then is an elliptic differential equation and, from the remarks 2, 3 and 
the initial condition (1.11), should have solutions 

ajiu, ^,v) = k sn{X u + F{i) - F{ri)) cd{\ u + F{i) - F{i])), (3.3) 



where k, as the module of eUiptic function, and A are two arbitrary constants and -F(^) 
an arbitrary function with the constriction k X = 1. 

Substituting (3.3) into the second of (1.14a) and its counterpart as well as the first 
and second of (2.26) and using elliptic function identities we have 

cd — sn +2mi{'r])k cd sn + u^ — Ui = 0, (3.4a) 

cd — sn — 2mi{'r])k cd sn + Uq — U4^ = 0, (3.4b) 

U1U4 + u^uq — cd — sn =0, (3.4c) 

uiUq + W4U5 — 2 m^{rj)k cd sn = 0, (3.4d) 

where m^{rj) and 771-1(77) are arbitrary functions satisfying 

7715(7?)^ -mi{r]f = -p. 

In the formulas (3.4)s and in what follows we simply write sn^ cd instead of elliptic func- 
tions sn{\u + F{S,) — F{r])) and cd{Xu + F{^) — F{ri)) for brevity. Using the symmetric 
operation we also have 

cd — sn +2 m4{^)k cd sn + u^ — u^ = 0. 
Since 777,1 (^) + iTT-ii^,) = one obtains 

cd — sn — 2rai{S,)k cd sn + u^ — u^ = 0. (3.4e) 

From (3.4c), (3.4d) and (3.4a) one also obtains 

— {cd + sn )ui + {cd — sn +2mi{ri)k sn cd.)u4^ + 2m^{r])k sn cdu^ = Q. (3.5a) 
If we do the symmetric operation with respect to the counterpart of (3.5a) then 

— {cd"^ + sn'^)u4 + {cd'^ — sr? — 2mi{^)k sn cd)ui + 2m^{^)k sn cdu^ = 0. (3.56) 
Similarly, one has 

(s7i + 2 mi{ri)k sn cd — cd )u^ — {cd + sn )uq + 2m5{r])k sn cdu4 = 0, (3.6a) 

(s77 +2mi{S,)k cd sn — cd )uq + {cd + sn )u^ + 2m^{^)k cd snu4 = {). (3.6b) 

Solving the equations (3.5a), (3.5b), (3.6a) and (3.6b) with respect to {ui,U4,u^,U(,} we 

have 

04(7^,^, 7?) _ H4 aQ{u,C,v) _ Hq 

ai{u,(,v) Hi' a5{u,^,r)) H5' 
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where 

Hi = (m.5(,f) + 7715(77)) cd + (777,1(^)7715(7/) + 777,5(^)777,1(7/)) /cs77, 

H4 = (777,5(,^) + 7775(7/)) Cd- (7771(^)7775(7/) + 7775(^)7771(7/)) /CS77, 

(3.6) 

-f^5 = ("7-5(0 +m-5(7/))s71 + (7771(^7775(7?) -7775(^7771(7/)) feed, 

-H'e = ("1-5(0 +m-5(7/))s77- imi{S,)mrj{r]) - mrj{6.)mi{r])) k cd. 
Let Til = HiX, U4 = H4X and 715 = iif5y, Tie = HqY. From (3.4a) and (3.4e) one obtains 

{Hl-Hl)X'^ =4 (7775(0 + "75(7?)) (7775(77)7771(0 +"75(0"7i (7/)) ksncdX^ 

= 2(7771(0 + mi{ri)k sncd, 

(3.7) 
(F| -HDy"^ =4 (7775(0 + "75(77)) (7775(77)7771(0 -m5(0"ii(^))^s"cdy2 

= 2(7771(0 ~ mi{rf))k sncd. 
It is easy to show using 7775(0^ — "t-i(0^ — 1/^^ that 



.2 _ 771l(0 +"T-l(7/) 



y2 



2("^5(0 + "i5(^))("ii(0"^5(^) + m^{x)mi{ri)) 

1 + k'^{m^{$)m^{vi) - 777i(0"T-i(?7)) _ -1 + k"^ {m^{i)m^{vi) + 777i(0"?-i(?7)) 

2("T.5(0 +"^5(??))^ 2(777i(0"T.5(??) +"1-5(0"il(^))^^^ 

7711 (0 -mi{ri) 



2(r775(0 + m5{ri)){mi{^)m5{ri) - mr,{x)mi{vi)) 
_ 1 + fc^(m-5(Om-5(??) + "T-i(0"T-i(^)) _ -1 + k'^{m5{i)m^{'n) - "t-i(0"^i(^)) 

2(w.5(0 +"1-5(??))^ 2(7r7i(OW'5(??) -"1-5(0"T-l(^))^^^ 

Hence gauge solutions of the colored Yang-Baxter equation (1.1) should obey the following 
forms 

ai(T7,e,r/) = A{i,r,)cd{\u + F{0 - F{7j)) + B{C,v) sn{Xu + F(0 - F{7j)), 
a2{u,i,r]) =03(77,^,77) = 1, 

04(77, e,r/) = A{i,r,)cd{\u + F(0 - F{r,)) - B{C,v) sn{Xu + F(0 - ^(77)), 

(3.8) 
a5iu,C,v)=C{tv)sn{Xu + F{C)-F{7])) + D{tv)cd{Xu + F{0-F{v)), 

aeiu,C,7]) = C(e,7?)s77(A77 + F(0 - F{r])) - D{tv) cd{Xu + F{0 - ^(77)), 
07(77, e, V) = ±fc sr7(A u + F(0 - i^(?7)) cd(A n + F(0 - i^(r/)), 
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where A;, as the module of the ehiptic functions, is an arbitrary constant and 



Mtv) = V{^ + G{0G{v)-H{0H{v))/2, 
B{i,v) = V(-l + G(OG(7?)+F(eMr?))/2, 
C{i,v) = 5^{l + G{OG{v)+H{i)H{ri))/2, 
D{tv) = V(-l + G(0G(7?) - H{OH{7j))/2 M, 

where 6^ = 1 and 

^ _ HiOGirj) - GiOHiv) 



(3.9) 



and G{^), H(^) satisfy G(^)^ — H(^)'^ = 1. If we consider the solution transformation D 
then the restrictive condition k X = 1 can be cancelled, namely, A 7^ is also an arbitrary 
constant. 

When the module k = 1 the Jacobian elliptic functions cd and sn should degenerate 
into 1 and tanh. Hence, we have 

ai(n,e,r/) = A(e,r?) + B{i,i^)tanh{\u + F{i) - F(r/)), 

a2{u,i,'n) = as{u,(,r]) = 1, 

a4(n,e,r/) = Aitv) - B{i,i^)tanh{\u + F{i) - F(r/)), 

(3.10) 
a^{u,i,r,)=C{i,r,)tanh{\u + F{i)-F{r,))+D{i,r,), 

aGiu,C,v) = C{tv)tanh{Xu + F{0 - F{r])) - D{tv), 

ajiu, e, ri) = ±tanh{X u + F(^) - F(r/)), 

where A{^,ri), B{S,,ri), G{^,rj) and D{^,rj) are defined by (3.9) and A 7^ is an arbitrary 
constant. 

If m^{tff' — 7711(77)^ = then the differential equation (2.29) can be rewritten as 
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So, following the calculation of (3.8) we can show that the gauge solutions are 

"''"■^■''> =^ ( co.(A°'fF(0 -V(,,)) +^g«'g''''""<^'' + ^«' -^W' 

a2{u,(,'n) = as{u,^,r]) = 1, 

a,{u,(,v) = X ( /^^f+J;^^^ - 2G{0G{r^)szn{Xn + F{0 - F{r^))) , 

\cos{X u + F(^) - F(ri)) J 

»'("•«■•'> = '' ( co,(a°+'f"(0-F(,)) + 2G(00(,).in(A„ + Fg) - F(„))) , 

°«<"-«-"> = ^ (- co^IA^+ViO-'fI,)) + 2G(0G(,),.„(A„ + F«) - F(,rt)) , 

a7(u,f,.,) = ±tan(Au + F({) - F(r,)), 

(3.11) 

where 

X = , ^ Y = ±X (3.12) 

2VG(OG(7?)' ^ ' 

and G(.^) is an arbitrary function. 

(3.2b) For the sub-case of m^{S) = —mQ{^) 

In the case of m^{ri) = —mQ{r]) weight functions of gauge solutions are 

cosh{Xu + F{0-F{7])) 

ai{u,^,r]) = a4{u,^,ri) = — — , 

cos{fiu + G{t,) - G(ri)) 

a2{u,C,r]) = as{u,(,r]) = 1, 

^ sinh{Xu + F{0-F{r])) 

a5{u,i.,V) = -a6{u,S„r]) = ± -— , 

cos(;UU + G(4) - G{ri)) 

ariu, ^, 7]) = ±tan{n u + G{() - G{r])), 

where A and ^ are two arbitrary constants, but not zero simultaneously, and F{^), G{S,) 
are two arbitrary functions. 

To prove it we first consider the sub-case of a 7^ 0. Then it follows that 

ai{u,(,r])a6{u,(,r]) + aA{u,^,'q)a^{u,^,ri) = 

by (2.26) and 

ai(u,^,r/)a5(u,C,r/) + ai{u,i,'q)aQ{u,i,'q) = 

by the symmetric operation. Since ai(u,S,,rj) ^ —a4,iu,S,,ri) owing to the initial condition 
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(3.13) 



(1.11). So one can get 

ab{u,i,i) = -ae{u,S,,r]), (3.14) 

ai{u,^,7]f -a5{u,^,r]f = 1 + 07(^,^,7/)^ 

and then mi{rj) = because of the first of (3.14) and the condition mi(^) + w,4(^) = 0. 
Then (2.29) has solution 

ariu, C, v) = tan{fi u + G{() - G{ri)) 

and hence (3.13) is true. 

Second we consider the sub-case of a = 0. Then it follows from the argument of 
proposition 2.2 in section two that the weight functions are 

ai{u, ^, 7]) = Ai{^, rfjcoshu — Aq{^^ rfjsinhu, 

a4{u,^,ri) = A4{^,ri)coshu + A^{^,ri)sinhu, 

a5{u,^,ri) = Arj{^,ri)coshu + A4{(,,ri)sinhu, (3.15) 

aQ{u, ^, 7]) = Aq{^, r])coshu — Ai(^, rfjsinhu, 

where Ai{$^,rj) = ai(0, ^, r/)(i = 1,4,5,6,7) are some functions with respect to colored 
parameters $,,rj to be determinate. It is clear that, Ai($^,rj){i = 1,4,5,6,7) satisfy the 
pure the colored Yang-Baxter (1.2). 

If we substituting the initial condition (1.11) into the ones obtained by putting v = —u 
and A = r/ in the system of (1.10) and combining resulting equation with the free fermion 
condition (2.12), one can show 

ab{u, S,, v) = -a5{-u, 7], 0, aaiu, S,, r?) = -aQ{-u, r], ^), 

(3.16) 
a7{u,^,v) = -a7{-u,ri,0, 04(^,^,1?) = ai{-u,ri,^). 

It is easy by (3.16) to show that 

A,{C,7])=A,{C,7]), Aeitv) = -MC,V)- (3.17) 

As mentioned in remark 2, all formulas obtained in this section and first two sections 
should hold for the pure the colored Yang-Baxter equation (1.2) except those obtained by 

23 



using symmetric operation. So we still should have 

1 + Aj{i, r,f - Alii, r,)A^{i, r?) - A^{i, r,)A(,{i, r?) = 0, 

/7(??)(Ai(^,r?)^6(e,^) +^4(e,r/)A5(^,r?)) = (l^{i^) +l^{^))Aj{i,r^), (3.18) 

/7(7?)(Ai(e,r?)2+A(C,r?)2-^4(e,r/)2-A5(e,r/)2=4/i(7?M7(e,r/), 

d 

where, as mentioned in remark 2, /«(??), (i = 1,4,5,6,7) mean ^^i(Ci'^)|f=r,- Substituting 

(3.17) into (3.18) we see 

Un) + k{v) = ^, hir]) = 0. 

Therefore, as we did for getting (2.29) we also have 

{-^AriC, v)f = /7(r/)(l + 2 ajiu, ^ V? + ^riu, t vf) (3.19) 

and then the solution for Aj(^, ry) is 

cosh{F{i) - F{r])) 



Ai{i,ri) = Ai{i,ri) 



cos{G{0 - G(7?)) 



A5(e,.)=-^a(e,.)= ""yi;^?-5/t^ ^'-'"^ 

cos{G{Q - Git])) 
Aj{tv)=tan{G{0-Gir])), 

where F and G are two arbitrary functions of single variable. Substituting (3.20) into 
(3.15) one can say (3.13) is also true for the case of a = 0, only fi = 0. 

Description above tell us that for the case of Hamiltonian coefficients m^{rj) = —mQ{r]) 
weight functions of a gauge eight-vertex type solution must be (3.13). 

Finally, straightforward calculation and computer symbolic computation can verify the 
following theorem. 

THEOREM 3.1 : Gauge eight- vertex solutions of the colored Yang-Baxter equation 
(1.1) are composed of (3.1), (3.8) and (3.13) and their degenerate forms (3.2), (3.10) and 
(3.11). 

The solutions (3.1) and (3.2) are called Baxter type solutions. They are just the 
solutions for "zero field" eight-vertex model by Baxter [3]. The solutions (3.8), (3.13) 
and their degenerate forms (3.10), (3.11) satisfy free fermion condition and are called 
Free-Fermion type solutions. 
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If we take A = 1, G{S,) = cosh{2S,), H{^) = sinh{2C) and F{^) = in solution (3.8) 
then (3.8) will reduce to the following solution 

ai{u,^,7]) = cosh{^ — ri)cd{u) + sinh{^ + r])sn{u), 

a4{u, ^, rj) = cosh{^ — r])cd{u) — sinh{(, + r])sn{u), 

0'5{u, ^, 7]) = cosh{^ + ri)sn{u) — sinh{^ — ri)cd{u), 

O'siu, S,, rj) = cosh{^ + ri)sn{u) + sinh{^ — ri)cd{u), 

ai{u, S,,T]) = k sn{u) cd{u) 
which is given in Ref. 20. 

§4. General solutions 

In this paper we have shown and classified all gauge eight-vertex solutions of the 
colored Yang-Baxter equation (1.1). These gauge solutions and trivial solutions (1.12a), 
(1.12b) together with five solution transformations discussed in the first section will give 
all eight-vertex type solutions. 

If we take in (3.8) and (3.9) 

and the solution transformation B with 

9{u,(,v) = Je{C)e{r])sn{^)sn{r]) , 

* sn[u/2) 

where the elliptic exponential 

e(C) = cn(C) +isn{C,) 
then using addition theorems for elliptic functions sn{C,) , cn{C,) and dn{() we can obtain 
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the following solution given in Ref. 12 

aiiu,^,-!]) = l-e(n)e(C)e(7?), 

(1 _ e(n)) 

sn{u/2) 

ai{u,e,,r]) =e{u) -e{S,)e{r]), 
a5{u,^,v) =e{0 -e{u)e{r]), 
a(i{u,i,r]) =e{ri) -e(n)e(C), 

aiiu, C, v) = as{u, (, rj) = -ik y/e{()e{r])sn{^)sn{r]) (1 - e{u)) sn{u/2), 

the detail of calculation of which is omitted. 

Similarly, all non-trivial general solutions can be also classified into two types. The 
first are Baxter type solutions if they can be obtained via gauge Baxter solutions and some 
solution transformations. The second are Free-Fermion solutions if they can be obtained 
via gauge Free-Fermion solutions and some solution transformations. 

According to the standard method by Baxter, for a given R-matrix the spin-chain 
Hamiltonian is generally of the following form: 

N ^ 

H = Y.^Jx CTjaj^l + Jy (j](j]+i + Jz CT|a|+i + -h (a| + a|+i)), 
i = l ^ 

where o"^, cr^, a'^ are Pauli matrices and the coupling constants are 

1 1 

Jx = jirn^ + mQ + mj + ms), Jy = -{m^ + me - nij - ms), 

Jz = t("T'1 - ms + 1714 - 7712), h = j{mi -1713-1714 + 711,2). 

In this paper we have proved that the hamiltonian coefficients of a gauge solution must 

obey 

22 22 

W-l = "1,4, 771,5 — "^6- 

It follows from the solution transformations B and D that 

(7711 — 7713) = (77I4 — 7712) , 7775 = TTlg 

for general solutions. This clearly describes the relation between classifications of eight- 
vertex type solutions and spin-chain Hamiltonians. For example, if Jx + Jy = h, J^ = 0, i.e. 
a special free-fermion model in a magnetic field, then one has 777,5 = "t-i ""t-s = —7714 + 777,2. 
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The corresponding solution of the colored Yang-Baxter equation should be (3.11). Then 
the transfer matrix is of trigonometric function type. 
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